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Abstract 

The simplest gravity duals for quantum critical theories with z — 2 'Lifshitz' scale 
invariance admit a marginally relevant deformation. Generic black holes in the bulk 
describe the field theory with a dynamically generated momentum scale A as well as 
finite temperature T. We describe the thermodynamics of these black holes in the 
quantum critical regime where T>A 2 . The deformation changes the asymptotics of 
the spacctime mildly and leads to intricate UV sensitivities of the theory which we 
control perturbatively in A 2 /T. 



1 Introduction 



At distance scales large compared to the lattice spacing, condensed matter systems are 
often described by quantum field theories (e.g. [1]). When the underlying electronic ground 
state undergoes a continuous nonanalytic change as a function of an external parameter, 
such as pressure, doping or magnetic field, the resulting quantum field theory becomes 
quantum critical [2] . Quantum critical theories are of interest both because they can control 
nonstandard phenomenology over large regions of the phase diagram and also because they 
provide an important theoretical starting point that is not necessarily weakly coupled. 

The quantum critical theories arising in condensed matter systems are scale invariant, 
but in general space and time need not scale equally [3 . The dynamical critical exponent 
z determines the relative scaling in which t — ► X z t but x — > Ax. Given recent excitement 
about application of the holographic correspondence [3] to condensed matter systems in 
general and quantum critical systems in particular [3 El El IE] j it was natural to ask whether 
the correspondence could be extended to 'nonrelativistic', i.e. 2^1, scale invariance. 

Following the logic of the original holographic correspondence, the authors of [§] wrote 
down a spacetime metric in which the scaling symmetry was realised geometrically: 



We will focus in this paper on a 3+1 dimensional bulk, with a 2+1 dimensional (putative) 
dual field theory. For z = 1 this spacetime becomes Anti-de Sitter space and the usual 
AdS/CFT correspondence is recovered. The full symmetry algebra of the spacetime ([I]) is 
often called the Lifshitz algebra [10] . Retarded Green's functions of operators dual to scalar 
fields in the bulk metric ([TJ were computed in (9] and found to have the form expected for 
field theories with Lifshitz symmetry. 

In order to obtain geometries like ([!]), with anisotropic scaling of space and time, it 
is clear that an anisotropic energy-momentum tensor is needed to source the gravitational 
field. The minimal way to achieve this is to include a vector field. It was shown in [TT] that 
a massive vector, i.e. Proca [12] , field is necessary. The Proca field is hodge dual to the two 
and three form action originally used in [9]. In this paper we will use an Einstein-Proca 
action. 

Given that additional structure beyond pure Einstein gravity is necessary to have a Lif- 
shitz metric as the ground state, it is important to understand the physics that is necessarily 
included in the theory with the addition of a Proca field. In the holographic correspon- 
dence, bulk fields correspond to operators in the dual field theory. The inclusion of a bulk 




(1) 
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Proca field indicates that Lifshitz field theories with gravitational duals come canonically 
equipped with a preferred vector operator, call it J a (not a conserved current). This is in 
addition to the energy-momentum tensor T ah which is dual to the bulk metric. 

The essential physics of the Proca field in the case z = 2 was explained in [9] : the Proca 
field (mixing with the graviton) describes a marginally relevant operator in the quantum 
critical theory. Turning on this operator induces a flow from the z = 2 theory to a rela- 
tivistic z = 1 infrared fixed point. The main objective of this present paper is to describe 
the renormalisation group flow more precisely and in addition explore the physics of the 
marginally relevant operator at finite temperature. That is, we are studying a canonical 
deformation away from criticality and looking at the interplay of the deformation scale and 
temperature scale. 

The value z = 2 is of particular interest both physically and mathematically. Physically, 
this value arises at 'deconfined' quantum critical points separating different Valence Bond 
Solid (VBS) phases, e.g. [131 HH EH HE]. The value z = 2 was also recently shown to 
arise in a simple string theory construction involving a distribution of e.g. DO branes 
|17j . Mathematically, z = 2 in the Einstein-Proca system is a special value at which the 
additional operator due to the Proca field becomes marginal. For z > 2 the operator 
is relevant while for z < 2 it is irrelevant |18[ [19]. In the marginal case, which in fact 
turns out to be marginally relevant, one finds a dynamically generated scale and nontrivial 
logarithmic running of various quantities. The renormalisation flow of the bulk fields has 
similarities to that encountered in 'improved holographic QCD' [20 j and in topologically 
massive gravity [21] . The holographic renormalisation of theories with marginally relevant 
operators is significantly more intricate than the relevant case that is usually studied; the 
asymptopia is mildly deformed away from its Lifshitz behaviour, leading to some interesting 
technical challenges which we partially resolve. 

In sections 2 to 4 below we introduce the Einstein-Proca theory and the space of solutions 
we will be considering. We characterise the asymptotic behaviour of these solutions and 
show how to compute finite thermodynamic quantities using holographic renormalisation. 
In section 5 we obtain numerical (planar) black hole solutions and use our previously devel- 
oped formalism to extract their thermodynamics. Our numerical results are summarised in 
graphs showing the dependence of thermodynamic variables on the ratio of the dynamically 
generated scale A over the temperature. 



2 



2 Einstein-Proca theory and Lifshitz spacetimes 



We will use the following action for a massive vector field coupled to gravity 
S 

with F = dA. The equations of motion following from this action are 







i 








~9* 


> 2 + ^ 2 



2k 



and 



c 2 



A u = 0. 



(2) 



(3) 



(4) 



The mass term has been tuned (relative to the cosmological constant term) in order for the 
Lifshitz metric with z = 2 to be a solution. With the normalisations chosen, the Lifshitz 
metric is 



ds 



Lif 



dt 2 dr 2 dx 2 + dy 2 
I 1 — — 

A 1 O O 



while the vector potential is 



A 



dt 



sf2 K r 2 

Both the metric and vector potential are invariant under the 'Lifshitz' scaling {t, x, y, r} 
{X 2 t, Xx, Xy, Xr}. With vanishing vector potential we also have the AdS^ solution 

dt 2 dr 2 dx 2 + dy 2 



(5) 
(6) 



ds 



AdS 



(7) 



It is important to note that there is no gauge invariance of the action ([2]). Therefore, the 
dual Lifshitz field theory does not automatically include a global symmetry, and the Lifshitz 
background should not be thought of as 'charged'. Rather, one can think of the divergence 
of the Proca field ^ near the 'boundary' r = as setting the preferred frame of the dual 
field theory. 

In this paper we will be interested in renormalisation group flows and in finite temper- 
ature solutions of the Einstein-Proca theory. These involve breaking the scaling symmetry 
at low energies and therefore we can introduce explicit functions of r in the metric. Thus 
we take as our Ansatz 



ds 1 
A 



f(r)dt 2 + ^+p(r)(dx 2 + dy 2 



■ g h{r)dt . 



(8) 
0) 
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With this parametrisation, the Lifshitz solution has 



Lifshitz: f = ^, P=^, k= ^^> (10) 
while the AdS^ solution has 

AdS: / = p=^r~Vi , h = 0. (11) 

With the Ansatz (|8]), the Einstein-Proca equations of motion are equivalent to the 
following three nonlinear ODEs for {f(f),p(r),h(r)} 

4h 2 rp' ^ r 2 f'p' r 2 p' 2 r 2 p" Q 

/ P VP 2p 2 P 
12h 2 rf r 2 f 2 2r 2 /V r 2 p' 2 r 2 f" . . 

20 + — " T + IF "~ ""^ ~ T - (12) 

Ah 2 r 2 h' 2 r 2 fp' r 2 p' 2 

10 A — — = 

f f fp 2p 2 

Before attempting to solve these equations, it is convenient to make the following change 

of variables 

p(r) = eJ T 'aM'" 1 * , f( r ) = e^ m(s)s ~ lds , h(r) = k(r)y/f(r) . (13) 

This transformation has the virtue of isolating the scaling ambiguity in p and / due to 
the possibility of rescaling the coordinates {t, x, y}. We will postpone fixing this ambiguity 
until our numerical section 5 below. Furthermore, now only first derivatives appear in the 
equations of motion. A final simplification is achieved by introducing 

x = y^O + 16/c 2 - 4mq - 2q 2 , 

which leads to the equations of motion 

rx = -8k - qx , (14) 
rq' = 5-2A; 2 - 3 | 2 -^, (15) 

rk' = _^_^! + ^_* + ^!. (16) 
q q 4 2 8q 

In deriving the above equation we have made an assumption about a sign which states that 

j^At < 0, or equivalently the Proca field grows as the boundary is approached. Changing 

the choice of this sign maps solutions of the Einstein-Proca theory onto solutions with the 

opposite sign of Af. However our requirement that the sign of the derivative remains the 

same for all r is significant. A priori the sign could potentially change at some radius on a 

given solution. We will find that this does not occur on the solutions of interest. 



4 



The equations ( 14 ) to ( 16 ) are, unsurprisingly, similar to those given in [9] for the Hodge 
dual theory. These are the equations we shall work with. In terms of these variables, the 
Lifshitz solution is now simply 



while AdS has 



Lifshitz: q 



AdS: 



-2 , x = 2V2 



k 



1 

71 



x = , k = 0. 



(17) 



(18) 



From the above equations we can construct an RG-invariant quantity. The combination 
1 



K 



fp\ — q + m + kx 



J Pf 4J0 + 16k 2 - 6q 2 + Aqx - x 2 



8q 



(19) 



is easily seen to be r- independent. This will turn out to be extremely useful in relating 
various physical quantities in our solutions. 



3 Asymptotic behaviour and the marginally relevant mode 

To understand the dual field theory physics of the various bulk fields, one starts by expand- 
ing near the 'boundary' of the spacetime at r = 0. The form of a general solution to the 



three first-order equations ( 14 ) - (|16|) as r — * contains three constants {A, a, f3}. One finds 



k Lfi H 1 | (-3 + A)-51og(-log(Ar)) 

y/2 V log(Ar) 21og 2 (Ar) 

o , 2 , A-51og(-log(Ar)) \ 

V l + log(Ar) + log 2 (Ar) + "V 

!(A0 4 log 2 (A0 f/J ( 1 + -i + ^og( log(Ar) ) + \ , i + _ _ x \ + ^ 



log(Ar) y \log(Ar 

where the ellipses denote the terms suppressed by a factor of ^i^hct^^ or more. Notice 
that the above are UV expansions in that they are valid near the boundary when Ar <C 1. 
We have not linearised the equations. 

Some comments about the parameters in the above expressions are in order. It appears 
that, apart from the three parameters {A, a, /?} controlling the three modes, there is an extra 
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parameter A that enters the solution. In fact the presence of A does not signal another degree 
of freedom in our system. Rather it is related to a 'gauge choice' in defining the scale A. 
More precisely, the solutions are mapped on to each other by the following transformation 

k(Ar; a, /?; A) = k{e x ' /2 Ar; e~ 2X ' {a - A'/?), e~ 2X ' f3; A + A') , (21) 

with the same holding for the other two functions q, x. This fact can be seen by verifying 



that all the series in the above solutions (20) can be obtained by putting A to zero and 
at the same time replacing log(Ar) with log(e~ A//2 Ar) and then re-expanding the series 
under the assumption |log(Ar)| 3> |A|. In order to compare different solutions, we can fix 
this ambiguity by fixing the parameter A once and for all. Subsequently, the parameters 
{A, a, /3} have an unambiguous meaning and can be compared among different solutions. 
From now on we will take A = 0. 

We see that there is one mode with an inverse logarithmic falloff and then two modes 
that go to zero as r 4 times logarithms. If we count the overall scale for the metric mode gu, 



left unfixed in (13), then up to logarithms we can think of this asymptotic behaviour as two 
'non-normalisable' modes that tend to a constant together with two 'normalisable' modes 
that go to zero. In the holographic correspondence, each operator in the field theory will be 
dual to a pair of modes. The 'non-normalisable' mode specifies the source for the operator 
while the 'normalisable' mode determines the vacuum expectation value. The power of the 
normalisable falloff determines the scaling dimension A of the dual operator. 

The scaling dimension A = 4 implies that we have two marginal operators in z + 2 = 4 
effective spacetime dimensions of the field theory. One of these will be the energy density 
T a , dual to gu, which necessarily has mass dimension [T tt ] = z + 2 = 4. This mode is 
marginal for all z. The other operator J , however, is only marginal for z = 2 [18, [19]. The 



expansions above (20) indicate that there is a logarithmic mixing between the two marginal 
operators. 



The inverse logarithmic falloff in (20) is consistent with a marginally relevant operator, 
whose renormalisation group flow dynamically generates the energy scale A 2 (cf. [20J). 
Recall that A is a momentum, so A 2 is an energy in a theory with z = 2. The fact that 
the operator is marginally relevant means that the Lifshitz scaling will be approximately 
recovered in the far UV, at energies much greater than the dynamical energy scale A 2 . We 
will quantify this statement very shortly as it turns out to be a little subtle. Previous works 
on black holes in the Proca-Einstein theory have tuned the marginally relevant mode to zero, 
e.g. [ITJ [181 032 1211 [231 IM] • This is the choice to work precisely at the quantum critical point. 
At the critical point the holographic renormalisation of the theory is significantly easier, as 
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the slow falloff of the inverse logarithm term introduces many new volume divergences of 
the on-shell action. 

In this paper we wish to consider the more generic solution in which the marginally 
relevant coupling is allowed to run. This will ultimately allow the study of crossover scaling 
physics as a function of T/A 2 (cf. [2]). The complication is that it will be challenging 
to holographically renormalise the logarithmic running. We shall approach this problem 
by working perturbatively in an expansion in inverse powers of log -rj . Concretely, the 
temperature will be much larger than the dynamically generated scale. 

Before turning to the issue of holographic renormalisation, we should look more closely 



at the asymptopia of the spacetime. From the expansions (|20j) and (13) we get 

1 2 1 A , -10-101og(-log(Ar)) 
l9tt{p) = F Q 7— rj— 47— 7 1 + iTZ773 +' 



t 2 ™^' U (Ar) 4 log 4 (Ar) V ^g(Ar 

1 n (A - p2 log 2 (Ar) / 4 + 51og(-log(Ar)) \ 
p9M ~ ^0-^52-^+ bi(AO + "J ' ( ] 

where Fq,Pq are constants. In this expression we see that in our Fefferman-Graham-like 
coordinates the asymptotic metric has logarithmic deformations away from Lifshitz be- 
haviour at leading order [23] . The curvature invariants and spacetime volume element y/— g 
in contrast are found to tend to their Lifshitz values plus inverse logarithmic corrections. 
This raises the question of whether these solutions can be included in an 'asymptotically 
Lifshitz' holographic correspondence. This concept has not been fully defined yet (some 
initial developments appeared in {25]) and will be an important question for future work. 
Analogous, though weaker, slowly falling off modes due to marginal deformations in the 
better characterised asymptotically AdS setup can be found for instance in [26, 27 . In 
the following we will present various arguments suggesting that one can make sense of the 



asymptopia (22 ), such as the computation of finite expectation values for operators includ- 
ing the energy. Our underlying reason for believing this behaviour should be allowed is 



that the UV deformation ( 22 ) was found in |9j to lead to a renormalisation group flow into 
the IR with the properties one should expect of a flow generated by a (marginally) relevant 
operator. We will find below that the finite temperature physics is also strongly consistent 
with this picture; for instance the deformation becomes increasingly less important at high 
temperatures. Marginally relevant deformations do not render the UV theory ill-defined. 



The logarithmic scaling in (22) would therefore be a manifestation of the fact that at any 
arbitrarily large fixed energy there is still a logarithmic deviation from scale invariance, as 
occurs for instance in QCD. 

To disentangle various effects, and simultaneously connect with previous results at A = 
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0, suppose we are interested in the behaviour of the series expansions at some high energy 
scale £ 3> A. Following the usual renormalisation group logic we can introduce a sliding 
scale fi S> A such that fir < 1. Then, we have 



By expanding in 



log Ar 



1 



A* 

A S Mr ' 



log fir 



< 1 



(23) 



(24) 



we can read off the leading order running of the couplings of the theory under the renor- 
malisation group. What we are doing here is formally showing how to take the A — > limit, 
which does not commute with the r —* limit. We first expand in j^w keeping fir fixed. 
We can then take fir small if we wish. When we later turn to high temperatures, T A 2 , 
we will find that the re-expansion with T ~ fi 2 is accurate over a significant range of r. In 
order for this re-expansion to make sense, we need to keep the right quantities fixed as we 
turn off the dynamical scale A. Specifically we let 



a 







A* 
A* 



-log 



a+log (x) Mx)-^ oglog (x) + 



1 

12 







a — 2 log 



(25) 



This redefinition allows us to take consistently the fi/A — > oo limit of the expansions (20) 
with a, (3 and fir held fixed. 

This re-expansion, at the leading order, leads to exactly the linearized approximation 



to the differential equations (14)-(16) about the Lifshitz backround in the absense of the 
marginally relevant deformation, as was employed in earlier works such as |9l [T9j [T5] . Con- 



cretely, re-expanding our general solutions (20) in the way described above in terms of the 
new parameters leads to 

1 



k 



V2 
-2 - 

2^2 



+ (fir) 4 ((3log(fir) +«)+••• + 0{{firf) , 
- ^ r ) 4 (^(H + 5 -£)+••• + 0{{firf) , 



(fir) 4 01og(fir) + «)+••• + O^firf 



log logQVA) 



(26) 



The 



where the ellipsis denotes terms that are suppressed by factors of ~Xg~^/A)~' or iog(/x/A) • 
higher order terms can be found systematically. 

We can now return to the asymptotics of the geometry. After re-expanding the metric 



(22 ) in the above way, we can also rescale the t, x coordinates as 



after which the metric takes the form 

h-v - ?( i - s w + -)- (29) 

The Proca field behaves similarly, allowing us to conclude that in the regime where the 
re-expansion is valid the spacetime can be treated as asymptotically Lifshitz with paramet- 
rically small deviations. Because \i is arbitrary, we can perform the re-expansion at arbi- 
trarily high energies, consistent with the notion that at any fixed high energy the spacetime 
is Lifshitz plus a small breaking of scaling symmetry due to a marginally relevant operator. 
When we heat up the system, it will be natural to consider this argument with /i 2 ~ T. In 
the limit T/A 2 — > oo, precise Lifshitz invariance is recovered at energy scales much greater 
than the temperature. The various scales of interest are shown in figure [T] below. 



IR 



A 



uv 



is r 

A 



-1 



Figure 1: A few scales of significance in our discussion and solutions. The first one is A, which is 
the scale dynamically generated by the marginally relevant mode. In particular A = when the 
marginally relevant mode is turned off. For the high temperature black hole solutions we are able 
to discuss most explicitly, the scale A is behind the event horizon located at r7 . To the right is the 
sliding UV cut-off beyond which the strict Lifshitz asymptopia is destroyed (for a given choice of //, 
which can be taken arbitrarily large). 

One reason that the re-expansion is technically useful is that it moves the logarithms of r 



from the denominator in (20) to the numerator in (26). This is again familiar from QCD. To 
any fixed order in perturbation theory, logarithms of the scale in the renormalisation flow of 
the coupling will only appear in the numerator. However, because we know that the QCD 
coupling is asymptotically free, these logarithms must be resummed into the denominator 



so that the coupling flows to zero in the UV. Our original expansions (20 ) are the resummed 



expressions, while (26) has reintroduced a perturbation theory in the marginally relevant 



coupling by working at a scale \i which is much higher than the dynamically generated scale 
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A. Similarly, at temperatures T > A 2 we will find that the deformation can be treated 
perturbatively. 

Before moving on to discuss holographic renormalisation we briefly recall the physical 
role of relevant deformations in a weakly coupled realisation of z = 2 Lifshitz criticality and 
a difference with the present strongly coupled case. The critical Lifshitz action considered 
in |14| [15] [T6] is, up to 'instanton' terms and a marginally irrelevant coupling, 

S U f = \ J drd 2 x ((d r c/>) 2 + if (V 2 c/>) 2 ) . (30) 

This theory clearly has a z = 2 scale invariance. The above action admits a relevant 
deformation by the term piVcf)) 2 . If p is positive then the theory flows to a Lorentz invariant 
theory in the IR. This is very analogous to the renormalisation group flow found numerically 
in j9] and that we are investigating in more depth in this paper. The difference between 
the two theories is that in the strongly coupled theory with a holographic dual, the natural 
deformation is marginally relevant rather than simply relevant as in the weakly coupled 
case. If p is negative then one obtains a 'tilted' phase for the 'height field' 0, i.e. an 
expectation value for V</>, that spontaneously breaks spatial isotropy. The Lifshitz theory 



(30) is then understood as the quantum critical theory separating tilted and untilted phases 
at zero temperature. In this paper we are only investigating the approach to criticality from 
one phase. We hope our results will provide a framework for a holographic realisation of 
the full quantum phase transition. 

While we expect our bulk theory to be dual to a field theory with the same scaling 



symmetry as the above theory (30), we do not expect it to be dual to a theory with such 
a simple Lagrangian with a single scalar field. In the present paper we will not attempt to 
construct the dual field theory. Rather it will suffice to point out that there is a natural 
family of theories that has the counter-part of a matrix large N limit. Generalisations 
of this theory are candidates for the field theory dual to the bulk theory discussed in the 
present paper. To motivate the existence of such a large N theory, it will be illuminating to 



rewrite the Lagrangian (30) in the following way. Recall that a scalar is dual to a one- form 
in (2+l)-dimensions. Writing 

*d(j) = dA , 
the above Lagrangian is transformed into 

Sim =\ j drd 2 x (|B| 2 + K\V x E| 2 ) . 

The SU(N), SO(N), Sp(N) generalisations of the above U(l) Lagrangian have been con- 
structed in [28] and it would be interesting to establish the possible duality between such 
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theories and the various bulk theories, e.g. [91 |29], which produce the Lifshitz spacetime 
^ as a solution. 



4 Holographic renormalisation 

Using the equations of motion (12) it is possible to show that the bulk action ^ evaluated 
on-shell can be written as the integral of a total derivative and hence becomes a boundary 
term. Specifically, the Euclidean on-shell action gives the semiclassical free energy density 

K^'iSS^- (31) 

By density we mean the density per unit boundary volume. That is, the free energy is 
F = jd 2 x T, with no integral over the Euclidean time circle. 

Before evaluating this action, we will need to add the boundary action and countert- 
erms. Firstly there is the usual Gibbons-Hawking term required for the desired variational 
principle, which in Euclidean space evaluates to 

Jg_ h = \\imJ^K= (-) lim r(v7p)' , (32) 

where j a b is the induced metric on the boundary and n^dn is an outward-pointing unit 
vector to the boundary which defines the extrinsic curvature K^ u = — V^n^. The bulk 
and Gibbons-Hawking terms do not lead to a finite on-shell action, and therefore boundary 
counterterms are required to render the theory well defined. 

There has been some confusion in the literature over the correct holographic renormal- 
isation of asymptotically Lifshitz solutions to Einstein-Proca theory. The physics seems 
clear: the solutions describe two deformations of the Lifshitz theory by the operators T tt 
and J 1 and their corresponding expectation values. As neither of these operators appear 
to be irrelevant, they do not destroy the UV definition of the theory. Therefore we should 
be able to consistently describe the physics of all of these modes. It becomes a question of 
finding the correct boundary counterterms. 

The approach we take is similar to [29]. In some regards our situation will be easier 
than the situation in that paper, because the theory studied there leads to unusually slowly 
falling off modes^ even without the marginally relevant operator. Just as in that work, 
in order to preserve boundary covariance, we would like to consider counterterms in the 
form of a power series in A 2 = A^A^. See [19] for an alternative choice of counterterms 



1 Specifically, [29j have metric modes falling off like r 2 as opposed to r 4 in the coordinates |HJ). 
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involving the combination W 'A^A^, which does not seem to generalise straightforwardly to 
include the marginally relevant deformation. As we will explain in more detail shortly, it 
turns out that just three such terms, ^/7(co + c\A 2 + C2A 4 ), are enough to render all the 
physical quantities of interest to us finite. The complication due to the marginally relevant 
deformation is that the coefficients {q} become functions of log Ar. 

Consider the following boundary counterterms in the free energy density 




We chose to write the series in this form for convenience, since —^rA 2 — h = when 
evaluated on the Lifshitz solution. We will shortly write down the coefficients {cj} that 
make the total semiclassical free energy density 

.F = ^ E + ^G-H + ^c.t. = ^£) 2 lim^ , (34) 

finite. 

Besides rendering the on-shell action finite, the counterterms should also lead to a well- 
defined variational problem and in particular finite on-shell expectation values and charges 
obtained by variation of the action with respect to boundary fields. Naively, following the 
usual procedure we would write the on-shell variation as 

6F = ^-r ab 5 lab + J a 8A a , (35) 

but some care is needed in the present case. 

As shown in [30 in the more established AdS/CFT context, the usual charge defined 
by f d 2 x^/a^ a ki)T ab , where £, a d a is a boundary Killing vector, \fa = y/jxxlyy is the spatial 
volume element and k a d a is the unit normal to the boundary Cauchy surface, say t = const., 
is not conserved when non-scalar matter fields are present. Instead, following [I30"lll9|. when 
defining the energy-momentum tensor we should hold fixed the matter fields in the boundary 
vielbein frame, defined by the usual relation 

lab = V & l44 , 7/ = diag(±l,l,l) , (36) 

where we have used the indices a,b, ... to denote the local coordinates on the boundary 
and the hatted indices to denote the local tangent space coordinates. The boundary stress 
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tensor T ab is then given by 

bT = T\ be h a + J & 8A & , T ab = 7 a h e hh . (37) 
A short computation shows that 

T ab = T ab + J_ j(a A b) ^ 

>/7 



where the bracket denotes symmetrisation of indices. In particular, the energy density is 
given by 

£ = V° K & T ab = ^7 A + J' A t . (39) 

While the J 1 At term here is closely analogous to the standard chemical potential times 
charge term appearing in the canonical ensemble, we should remember that J 1 in our 
system is not a conserved charge. 
Taking the local frame to be 

e* = e l a dx a = s/fdr , e x = ^ppdx , e s = y/pdy , (40) 

the quantities of particular interest can be computed to give 



2 S - ' (A'"'' A'-."'' 

lab K 2 V 



r afe = — ^- = ^ ^ ab -A- 7 a6 (41) 

y/l^ab K 2 V / 



+ ^E-(n-^S) +2 V^(-7^-5) )• (42) 



and 



2 _ l\Tl-l 

n=0 



1* 

.9 



ilirnJfpl-lx + kJ^ncnik 2 -^- 1 ) , (43) 

K r— >0 \ Z ' — » / 

\ n=0 / 



and all other ^7° = 0. We have used the equation of motion (16) to obtain the second line 



of the previous equation. Putting these expressions together, we obtain 

^{^^(v-Y^t^-^) ' <44) 

It will be very useful to note that the above expression for the energy density is related to 
the free energy in a simple way. Observe that 



fp 



rp' xk | 1 
~p 2 



+ ^X>(^ 2 - hr) = \^fp{ r -T + r -^ + jz ^ 2 - ir) + K . ( 45 ) 

n=0 / \ J P n=Q J 
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where K is the RG-invariant quantity we found in (19). Evaluating both sides in the UV 



(r — ► 0) and using (44) and (34) leads to the simple relationship 

£ = ^)V (46, 

As we shall see later, this gives exactly the 'integrated first law of thermodynamics' when 
considering a black hole solution. Finally we can easily compute the pressure and find 
V = —y/jr x x = —T as expected on general thermodynamic grounds. The off-diagonal 
components such as T l x = T t x vanish on our solution. 

As we have implied in the preceding discussion, in order to renormalise the on-shell action 
we will need three boundary terms given by three non- vanishing coefficients co,ci,C2- It 
turns out that instead of being constant, they are series in log (\ r ) ; with r understood to be 
evaluated at a small cutoff value. This explicit dependence of the boundary counterterms 
on the cutoff is expected: the marginally relevant perturbation of the theory breaks scaling 
invariance at any fixed scale. This breaking of the scaling symmetry can be viewed as a 
conformal anomaly [31]. However, scale invariance should be recovered in the strict r — > 
limit, or equivalently the A — > limit, and hence we can expect inverse powers of log Ar to 
appear. Imposing finiteness of the free energy and energy densities we find 



65 + 1201og(-log(rA)) 
31og^(rA) ' 9 log 3 (r A) 



-91-1201og(-log(rA)) 
1 + 31og(rA) + ~ 91og 2 (rA) { ' 

5 35 a- ^log(-log(rA)) 
° 2 ~ ~3 + 91og(rA) + log 2 (rA) ' 

where the ellipses again denote terms suppressed by extra factors of log(— log Ar)/log Ar 
or 1/logAr. A few remarks about these coefficients are in order. Firstly, notice that the 
first two coefficients are given in terms of infinite series, while the series expansion of the 
third coefficient terminates at l/log 2 (Ar). This can be understood in the following way. 
Observe that the divergences of the free energy and energy before renormalisation are both 
given in terms of r™ 4 multiplied by an infinite series in 1/logAr, plus a finite number of 
divergent terms of the form log n (Ar) for some positive integer n. Hence we have altogether 
two infinite series and a few extra terms which need to be cancelled by the counterterms. 
Unsurprisingly, they are exactly cancelled by counterterms given by two infinite series and 
three more terms. Relatedly, there is an ambiguity in the final counterterm given by a 
real number denoted by a in the above equation. This ambiguity does not afflict any of 
the thermodynamical quantities we shall study later. Nevertheless it does represent an 
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ambiguity in the quantity J^, which reflects the extent to which the system reacts upon 
changing the boundary (i.e. background) value of the Proca field At. Such ambiguities 
are physical and common in the presence of logarithmically running couplings. It is also 
possible to add more subleading terms in the series for C2, but these terms will not affect 
our results for £ or J 1 . Finally, it is also possible to add higher power counterterms such 
as cs(k 2 — ^) 3 . Nevertheless we shall make the above minimal choice which will be sufficient 
for our purposes. 



After incorporating the appropriate counterterms given by (33) and (47), we obtain the 



following finite expressions for the physical quantities in terms of the parameters a, (3 in our 

2 V2 



expansions (20): 



£ 

J 1 



9 v 

K 

1 1 /_27 



-5/3 + 6a) , 
(7/3 + 6a ) , 

V 

+ 2a ) (3 . 



(48) 



From here we see that evaluated on the pure Lifshitz solution we have 



£ = T = J l 



0. 



Furthermore, recall that we have an RG-invariant quantity (19) which can be written in 
terms of the parameters a, (5 as 



K 



1 



fp\ — q + rn + kx 



~9~ 



(£ + 6a) , 



(49) 



consistent with the relation (46) between the quantities £, T and K. 



5 Finite temperature 

Several works have studied black holes in Einstein-Proca theory for various values of z, 
e-g- [13 [22l ISSl HI]- As we noted above, these studies focused on the quantum critical 
theory and tuned the marginally relevant deformation to zero. We will begin similarly to 
those studies, by expanding the (planar) black hole solution near the horizon, and then 
numerically integrating outwards towards the boundary. Since we have now learned how to 
analyze solutions which contain the marginally relevant mode at the boundary, we will be 
able to explore the energetics of these black objects for a family of black holes with different 
fluxes. 
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5.1 Expansion and physical quantities near the horizon 



The black hole horizon is defined by f(r+) = 0. Given a horizon we can proceed to expand 
the solution near r = r + . In the coordinates we have chosen, regularity requires the gu 
component of our metric to have a double zero at the horizon, while the g xx component 



should go to a nonzero constant. Solving the equations of motion (12) in terms of our 
Ansatz d8|) under these requirements, we find the expansions 



f(r) 
p{r) 
h{r) 




Po 1 + 



2hl 



+ 



+ 
5 



11 + 16/i| 
~12 



2hl 



1 - 
r 



+ • 



+ • 



h n 1 



+ ho 1 



+• 



(50) 
(51) 
(52) 



Varying the constant ho produces the family of black holes we will study. Varying r+ 
only causes various quantities to scale in a way determined by their dimension, and so r + 
will drop out of the dimensionless quantites that we will compute. The constants /o and 
Po are unfixed at the horizon and determine the clock and ruler of the system. We will 
eventually fix them by imposing Lifshitz asymptotics on a scale set by the temperature 



T » A 2 as in (28) and (29). 



We can calculate a few physical quantities characterising the horizon of these black 
holes, specifically their temperature, entropy density, and horizon flux density. These are 
given by 



r=r+ 



r+ /ldV 
2vr V 2 dr 2 

2^(J p(r + ), 

£gr + ( p dh\ 

77^) 



r=rj r 



(53) 
(54) 
(55) 



Using the expansion (50) near the horizon, we obtain 

2 



Ik 



o 



2vr 



2irpo 



2h p 



(56) 



Additionally, we can compute the value of the r-independent quantity K from (19) 
above, which we can then rewrite in terms of T and s: 

'KN2 



K = \ffopo = Ts 



(57) 
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From (46) this immediately implies that 



T = £ -Ts, (58) 

which is the integrated form of the first law for these black holes. This result is an important 
consistency check for our framework and we will use it later as a check on the accuracy of 
our numerics. 

Before turning to numerics, we can make some analytic predictions about these physical 
quantities in the limit where the marginally relevant mode is not excited. In this case, when 
A = 0, we expect from Ward identities that the pressure should be equal to the energy, 
because our anisotropic scale invariance is unbroken. Recalling that T = — V , we therefore 



expect that J-q = —So- By examining the expressions for T and £ in (48), we see this 
relation can only be true if (5 = 0. Combining with the integrated first law, we thus expect 
the following results at A = 0: 

So = -To = \tso , Jl = 0. (59) 

This relation between £ and Ts at A = is a special case of the relationship for general 
values of z demonstrated in 



5.2 Integrating towards the Lifshitz Boundary 

Equipped as we now are with various physical expectations and mathematical relations, 
let us explore the results of numerical integration. Our basic technique is to take the near 



horizon expansions in (50 ) and then use the equations of motion in ( 12 ) to integrate towards 
the boundary. We will always use r/r+ as our variable, and we will want to explore unitless 
quantities such as T '/Ts, s/T and A 2 /T. 

When performing this integration, we find that we cannot use a value of ho bigger than 
/imax ~ .9714. At larger values of the flux, the numerical integration procedure produces 
metric functions which grow exponentially as we try to take r to 0, and we can never reach 
a boundary. This behaviour was first noted in |23j . While this implies a maximum value of 
the horizon flux for the black holes we consider, asymptotically we will see that this limiting 
value corresponds to the high temperature limit A 2 /T — > 0, in which the deformation is 
turned off and we recover a black hole in the pure Lifshitz spacetime. These are the black 
holes that have been considered in the previous works we mentioned above. 

We also have a minimum value of the flux; when the flux is strictly zero at the horizon, 
we reproduce standard asymptotically AdS black holes in which the flux is zero everywhere. 
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However, the limit in which the flux is taken to be very small is nontrivial. This limit 
corresponds to the zero temperature limit A 2 /T — > oo. If we take this limit with A fixed 
then we recover the RG flow from UV Lifshitz to IR AdS found in [9]. In units of the 
temperature, the crossover regime to Lifshitz scaling moves further and further away from 
the horizon. This is why if we keep the horizon quantity T fixed in this limit we simply 
recover pure AdS black holes. 

Thus we will find that tuning the horizon flux via ho corresponds to interpolating 
between the zero temperature RG flow of [9] and the black holes solutions without the 
marginally relevant deformation of e.g. jl H [TBI [T9| [22] 123) 124] . For the remainder of the 
paper, for technical reasons to be discussed shortly, we will be most interested in solutions 
whose flux is just below the maximum; that is, we will be exciting only a small amount of 
the marginally relevant mode. 



5.2.1 Matching A, /q and po 

Before we can attempt to obtain the thermodynamic variables of interest through formulae 



like (48), we need to characterise the asymptotics by extracting A, /q and p$. The first of 
these is a physical quantity in its own right, while the latter two will tell us how to rescale 
space and time in order to consistently compare between different solutions. 

The scale A can be found by fitting the numerical results to the asymptotic expectations. 



Concretely, we fit to the expansions in (20) sufficiently close to the boundary that the r 
terms are suppressed. For the solutions we consider, we find it is sufficient to evaluate near 
log r/r+ ~ — 10 5 . This procedure allows us to find the value of A for a given solution, 



provided we fix a value of A in (20). In practice the dimensionless quanitity we obtain is 
Ar + which we will be able to relate to A 2 /T once we have found /q. A quick check we 



performed on the numerics was to reproduce the 'covariance' behavior of (21) for A as a 
function of A. As in previous sections, we will now set A = 0. 

Upon extracting A in this manner as a function of Hq we find, that Hq — /i max gives A = 0. 
Therefore the maximum flux solution is also the solution with none of the marginally relevant 
mode excited. We now have a first order understanding of why problems arise when trying 
to integrate solutions with more flux; we would have to match them onto asymptotics with 
negative A. 



Extracting fo and po is more subtle. These divide out of the equations ( 12 ) and must be 
fixed by specifying the asymptotic normalisation of the metric components gu and g xx . The 



leading order logarithmic running to the boundary uncovered in (22 ) means that the r —* 
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limit of r A gtt and r 2 g xx remains r dependent. This was manifested in (27) where we needed 
to rescale time and space in a /^-dependent fashion in order to obtain the near-Lifshitz form 
of the metric (28) and (29). As we will be working close to A = 0, at small A 2 /T <C 1, we 
can take p ~ r^ 1 in (28) and (29). In this high temperature regime we thus have a well 
defined region in which the near-Lifshitz form holds. This will allow us to fix fo and po in 
the undeformed theory and then extend away from A = 0. 

In figure [2] we illustrate how the numerically generated / and p can be fitted over a finite 



range, at high temperatures, to the near-Lifshitz form (28) and (29). In fact we used these 
expansions to fifth order in l/log(Ar+). This allows us to obtain fo and po- The figure 
depicts the fitting procedure for ho = .962, which corresponds to logA 2 /T ~ —200. The 
fitting is even more robust for the smaller values of A which we will usually consider. 




Figure 2: Left: fr^/for 4 ^ versus logr/r + . The dashed red line shows the function fr^/for 4 as 
given in (28l, including higher order terms, for /qt+ = 27.72. Right: pr 2 /por^_ versus logr/r + , with 



the dotted line corresponding to po^+ = 1-89. The black line is the numerical curve and the black 
dots are the values used for fitting. 

Having found A, fo and po for a solution with a given ho, we can plot the (dimensionless) 
entropy density over temperature as a function of log(A 2 /T). This is shown in figure [3] In 
particular at A = we find 

s 4vr 2 /A 2 . . 

tL-tiU) • (60) 

We have only found this value numerically, to several digits of precision, and do not have 
an a priori understanding of its origin, since we do not know the form of the analytic 
solution. The value we find seems to agree with the value s/T pa 11 A quoted in |23 j if we 
set {£/k) 2 = \. 
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Figure 3: Entropy density over T versus logA 2 /T for ho — 0.9713 to ho = 0.9708. Dots are data 
points, which are joined by straight lines. 

5.2.2 Energy, Free Energy, and J 1 

Our next objective is to compute the temperature dependence of F/Ts, £/Ts and J 1 jTs 



with the numerators given in (34), (44), and (43) respectively. As well as being unitless, 



these ratios are independent of the constants /o and po and are therefore not sensitive to 
any, possibly delicate in this logarithmically running theory, choice of normalisation of time 
and space. 

We cannot directly extract a and (3 from the numerics and then use the expressions 



(48) for the energy etc., because the r terms in the near boundary expansions (20) are 



exponentially small compared to the logarithmically decaying modes. Instead we will use 



the expressions ( 34 ) , ( 44 ) , and ( 43 ) together with the series expansions for the counterterm 



coefficients in (47). In practice we only know the counterterm expansion (47) up to some 



fixed (arbitrarily high in principle) order. The main fact that restricts us to the high 
temperature regime is that the truncated counterterm series is more accurate there. With 



A 2 /T <C 1 we can estimate that truncating the expansions (47) at order l/\og N (Ar) will 
be reliable for 



log 



< N ^ 

< T log 



log 



A^ 



(61) 



We see that increasing the number N of terms allows us to move closer the boundary, as does 
increasing the temperature. At sufficiently large N and temperature we can expect to obtain 
a constant intermediate region between the near horizon IR effects and the asymptotic 
divergences. From this region we can read off our thermodynamic variables. This process 
is illustrated in the plots of figure |4j for which we expanded the counterterms to ten orders. 
There is an additional complication in obtaining J 1 . For T and £ the counterterms not 
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Figure 4: Plots of J-/Ts, £/Ts and J 1 jTs as a function of log r/r + for /i = .971, corresponding 
to log(A 2 /T) w —3210. The near horizon region is dominated by IR effects while near the boundary 



divergences set in due to use of a truncated series of counterterms ( 47 1 . The quantities are well 



defined in the intermediate region, which can be made large by working with a sufficiently high 
order counterterm expansion. Black lines are numerical results while the red dashed line is a fit to 
the analytical renormalised expressions. 

only remove all divergences, but also all subleading terms in 1/ log (r/r+). Therefore the 
only significant error is due to truncating the counterterm expansion. This is not true for 
J ; while the counterterms remove the divergences, there are still inverse powers of log r/r + 
that vanish asymptotically but which are nontrivial in the regime in which we can read off 



the value of J 1 . The behaviour can be derived analytically from (43) and so we can read off 
i7* accurately by fitting the numerics to an expansion of several orders in l/log(r/r + ) in 
the intermediate regime. Finally, we should recall that the quantity J % is dependent upon 
ambiguities in the counterterms and is not a conserved charge. 



5.3 Exploring the dependence of £ and T on logA 2 /T 

We finally have all of the tools available to compute T '/Ts and £/Ts as a function of 
logA 2 /r. The results are shown in figure [5] for a range of high temperatures. 




£ ".5006 

Ts 
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(A 2 




Figure 5: Plots of T ' /Ts and £/Ts , and T jE , as functions of log(A/T 2 ). The range of h is from 
0.9712 to 0.9698, corresponding to logA 2 /T from about -12000 to -1150. The dots are numerical 



results, while the lines are the fits in equations (63 1, (64 1, and (651. 
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The first observation to make about figure [5] is that we reproduce the anticipated results 
(59) as A — ► (that is, as logA 2 /T — > — oo). We find numerically 

^- = -- A = I ^° = -1 (62) 
s T 2' s T 2' So { ' 

Beyond this limit, we find stable fits for each of these quantities as an expansion in inverse 
powers of log A 2 /T. For T we find 

Ts = -\ + ^]f + -- (63) 

The coefficent of 1 in the numerator of the second term here is not known analytically, it is 
obtained to some precision from our numerics. The fit (63) is also shown in figure [5] The 
simplicity of this term may suggest the existence of a simple analytic solution. 



We can now use the integrated first law (58) to predict the expansion for S : 

S _ 1 t 1 

Ts ~ 2 + logA 2 /T + "' ' ( ' 

which, as Figure ([5J shows, is strongly supported by the numerical results. Note that ^ is 
completely determined by , even beyond the leading order terms written down in (63 64 ) . 



Combining the expansions (63) and (64), we expect J-/S to fit 

which is also seen to hold within numerical precision for this range of log Is? /T. 

Thus the graphs in figure [5] provide a check for the integrated first law, or rather that our 
numerics satisfy the law. In order to see how precisely our numerics satisfy this requirement, 
we can examine figure [6j which plots —T/Ts + £ /Ts — 1. Since we expect this quantity to 
be zero, its size gives us an estimate of our numerical error. As the plot shows, our error 
is less than 1 part in 10 5 for the range considered. While the variation of our variables of 
interest in figure [5] are also small, the variation within the range of A 2 /T we consider is at 
least two orders of magnitude bigger than the error. Furthermore the error does not appear 
to be systematic. 

Lastly, we present a plot of J 1 /Ts against log A 2 /T, in figure [Tj As we can see from the 
figure, we find that J 1 = at A — > as expected. We have not been able to fit this curve 
to a simple function. 
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Figure 7: J*" /Ts for ho from 0.9712 to 0.9698. Dots are numerical data points joined by lines. 

6 Discussion 



Our main concrete result is perhaps equation (63) expressing the A 2 /T depend ence of 
the free energy of the theory upon approach to Lifshitz quantum criticality at A = 0. It is 
important to characterise the strongly coupled physics of (marginally) relevant deformations 
away from criticality. To the extent that such deformations are generic, they should be 
included in attempts at holographic model building for condensed matter systems along the 
lines of [IT]. Beyond the thermodynamics we have studied in this paper, it will be of interest 
to consider correlators of fields in the deformed background and describe their behaviour 
upon approach to criticality. 

More generally, however, we have tried to grapple with the technical challenge of renor- 
malising the leading order logarithmic running of the deformation of the theory generating 
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the scale A. We have worked perturbatively about the A = point (equivalently, at high 
temperatures relative to A) and successfully computed finite thermodynamic quantities. 
Obtaining these results required a certain amount of numerical acrobatics and it would 
clearly be advantageous to have a better formal understanding of the spacetime asymp- 
topia as well as a more general approach that would allow computations all the way to the 
complementary T = limit. 

While the deformation dynamically generating the scale A has many of the hallmarks of 
a marginally relevant deformation of a Lifshitz invariant theory (e.g. at T = the theory 
runs to a new IR fixed point while for A 2 <T one recovers scale invariant results), we have 
not shown rigorously that our spacetimes can be considered 'asymptotically Lifshitz'. In 



fact, the leading order logarithms in Fefferman-Graham-like coordinates (22) may indicate 
that this is not strictly the case. The need to rescale space and time together with the 
energy scale in (27) may suggest a weak generalisation of the notion of a fixed point(^j We 
hope that future work will elucidate this question. 

In weakly coupled realisations of z = 2 Lifshitz symmetry, reviewed briefly at the end of 
section 3 above, the relevant operator that drives the theory to a relativistic IR fixed point 
is also responsible for the existence of 'tilted' phases breaking spatial isotropy. A fascinating 
open question is to identify such phases, if they exist, in the gravitational dual. One is led 
to wonder if such phases exist beyond the maximal horizon flux that we found. 

There are various further natural extensions of our work: An analysis similar to ours 
should be possible, perhaps easier, for z > 2 where the deformation mode becomes strictly 
relevant. It would be very helpful to have analytic solutions, we found some encouraging 
signs that this might be possible. Finally, logarithmic modes can occasionally make the 
stability of the spacetime quite a subtle question - it would be interesting to characterise 
the fluctuations about our marginally deformed backgrounds. 
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